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ABSTRACT 

Classical Cepheid variable stars are crucial calibrators of the cosmic distance 
scale thanks to a relation between their pulsation periods and luminosities. Their 
archetype, 6 Cephei, is an important calibrator for this relation. In this paper, 
we show that S Cephei is a spectroscopic binary based on newly-obtained high- 
precision radial velocities. We combine these new data with literature data to 
determine the orbit, which has period 2201 days, semi-amplitude 1.5kms“^, and 
high eccentricity (e = 0.647). We re-analyze Hippareos intermediate astrometric 
data to measure J Cephei’s parallax {w = 4.09 ± 0.16 mas) and hnd tentative 
evidence for an orbital signature, although we cannot claim detection. We esti¬ 
mate that Gaia will fully determine the astrometric orbit. Using the available 
information from spectroscopy, velocimetry, astrometry, and Geneva stellar evo¬ 
lution models (M^cep ~ 5.0 — 5.25 Mq), we constrain the companion mass to 
within 0.2 < M 2 < 1.2 Mq. We discuss the potential of ongoing and previ¬ 
ous interactions between the companion and S Cephei near pericenter passage, 
informing reported observations of circumstellar material and bow-shock. The 
orbit may have undergone signihcant changes due to a Kozai-Lidov mechanism 


^Swiss National Science Foundation Fellow 



2 


driven by the outer (visual and astrometric) companion HD 213307. Our discov¬ 
ery of 6 Cephei’s nature as a spectroscopic binary exposes a hidden companion 
and reveals a rich and dynamical history of the archetype of classical Cepheid 
variables. 

Subject headings: binaries: general, binaries: spectroscopic, stars: distances, 
stars: individual: 5 Cephei = HD 213306 = HIP 110991, stars: oscillations, 
stars: variables: Cepheids 


Introduction 


Classical Cepheid variable stars (from hereon: Cepheids) are precise Galactic and extra- 
galactic distance tracers and thus of crucial importance for cosmology. The prototype of this 
class of stars, 5 Cephei (HD 213306, HIP 110991), has been extensiveljj^ studied ever since 
the discovery of its variability 230 years ago by Goodricl^ ( |1786[ ). Until Baade’s proposed 
method ( Baade||1926 ) to test Shapley’s pulsation hypothesis ( Shapley||1914 ) bore fruit in the 
early to mid 20th century, Cepheids were thought to be binary stars on eccentric orbits (see 


e.g. Gautschy 1997), inspiring a great deal of research, not least that by Ghristian Doppler. 


The hrst scientists to measure the variability of 5 Gephei’s spectral lines and to determine 


their velocities were Belopolsky (1894, 1895) and Moore (1913), before radial velocities (RVs) 


became available for many Gepheids thanks to the observations by Joy (1937). Two decades 


later, Shane (1958) conducted a detailed analysis of 5 Gephei’s RV curve and concluded that 


no evidence of long-period changes in the mean velocity could be seen, a result that was 


conhrmed much later by Kovacs et al. (1990). 


are 


Nowadays, Gepheids are known to be pulsating variable stars, although many Gepheids 
also known to be binarie^ The binary fraction of Gepheids is being studied intensively 


with most recent estimates of the total binary fraction ranging around 60%, see e.g. Evans 
et al. (2013) and Szabados et al. (2013). However, Gepheids cannot reside in very close- 


in binary systems (e.g. Neilson et al. 2014) due to their nature as evolved (super-)giant 
stars, which results in observed minimum orbital periods of approximately one year. For 
long periods (> 10years), practical constraints create signihcant observational bias against 
companion detection. 


1 
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TMBAD (http;//simbad.u-strasbg.fr) lists more than 600 articles related to <5 Cephei. 
;f. the Cepheid binary database by 


Szabados 

(2003 

1 at 

http://www.konkoly.hu/CEP/nagytab3.html 
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To our knowledge, 6 Cephei has not been shown to be a spectroscopic binary prior to 
this work. 5 Cephei is, however, a known visual binary (Fernie 1966), whose companion 


HD 213307 (= HIP 110988) is itself is an astrometric binary and thought to be physically 
associated (see Benedict et ah[2002 and references therein). HD 213307 was proposed to also 
be a spectroscopic binary (Herbig, priv. comm, mentioned in Fernie||1966). Finally, 5 Cephei 


is usually considered to be a member of the loose association CepheusOBO (de Zeeuw et ah 


1999 

van Leeuwen et al. 

2007 

Majaess et al. 

2012 


Precise trigonometric parallax of 6 Cephei has been measured by [Perryman fc ES~^ 


(1997) 

van Leeuwen et al. 

(2007 

), and 

van Leeuwen et al.| ( 

2007 

) using observations made 

by the Hipparcos space mission and by 

Benedict et ah] ([2002 

) using measurements obtained 


with FGS 3 on board the Hubble space telescope {HSTj. Notable previous distance estimates 


include those reported by Fernley et al. (1989), Gieren et ah (1993), Gatewood et al. (1993), 


Mourard et al. (1997). 


Merand et al. (2005) employed infrared long-baseline interferometry to study the Baade- 


Wesselink projection factor of 6 Cephei and were later ( Merand et al.||2006 ) able to show the 
presence of an extended circumstellar envelope. The presence of this circumstellar environ¬ 


ment was conhrmed independently (using different methodologies) by Marengo et al. (2010) 


and Matthews et al. (2012) 


In this paper, we present the discovery of the spectroscopic binary nature of 6 Cephei. 
This discovery is demonstrated using new observations that are presented in Sect.[^ Using 
these new data we reveal the presence of a hidden companion in Sect. After combining 


our new RVs with literature data in Sect. 3.2, we determine the orbital solution in Sect. 3.3 


We re-analyze the Hipparcos intermediate astrometric data from the van Leeuwen (2007) re¬ 


duction in Sects. 4.1.1 and 4.1.2 to measure parallax and to investigate whether the available 
astrometric measurements are sensitive to the motion due to binarity. In Sect. 4.2, we inves¬ 


tigate Gaia's expected sensitivity to the binary motion. We discuss how our discovery helps 
to better interpret other observations and begins to draw a complex picture of 5 Cephei’s 
rich and dynamical history in Sect. [^before concluding in Sect.[^ 


New Hermes Observations 


Line-of-sight (radial) velocities were measured from 136 observations taken between 
September 2011 and September 2014 using the fiber-fed high-resolution (R ~ 85 000) spec¬ 


trograph Hermes (Raskin et al. 2011) at the Flemish 1.2m Mercator telescope located at 
Roque de los Muchachos Observatory, La Palma, Canary Islands. We utilize the high- 
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resolution fiber (HRF) mode for all observations, since this is the most commonly used, 
i.e., best-understood, observing mode available for Hermes. The HRF mode offers optimal 
efficiency and the highest available spectral resolution. 


The reduction pipeline available for Hermes performs pre- and overscan bias correction, 
flathelding using Halogen lamps, and background modelization, as well as cosmic ray removal. 
ThAr lamps are used for the wavelength calibration. RVs are determined via the cross¬ 
correlation technique (Baranne et ah 1996 Pepe et al.|2002 ) using a numerical mask designed 
for solar-like stars (optimized for spectral type G2). 

These observations were started as part of a search for Cepheids belonging to open 
clusters ( Anderson et al.||2013 ) with the goal of quantifying the precision limit for a classical 
Cepheid using Hermes RVs and of having high-quality spectra available for further study. 
The average signal-to-noise ratio of the spectra is higher than 200 near 6000 A, with some 
spectra reaching up to 400. 

Table[^ shows a sample of our Hermes RV measurements that we are making publicly 
available at the CD^ together with our standard star observations. 


2.1. Zero-point, Stability, and Precision of Hermes RVs 


To determine the RV zero-point of Hermes, we observed eight RV standard stars listed 
Udry et ah ( 1999a||b |p see Tab. i These standard stars were chosen to cover the range 


m 

of spectral types that 5 Cephei exhibits during its pulsation cycles. We thus determine a 
mean systematic offset of 55ms“^ with respect to the ELODIE and CORAVEL zero-point. 
We estimate this zero-point offset to be accurate to approximately 10ms“^ and note that 
additional scatter in the difference between Hermes and literature RVs can exist for various 
reasons of astrophysical origin, including binarity and planetary companions. We therefore 
increase our error margin by adding 10ms“^ in quadrature to the Hermes RV uncertainties 
when correcting for zero-point offsets. 

We investigate the long-term stability and precision of Hermes RVs using the RV stan¬ 
dard star observations. Thanks to the very high signal-to-noise ratio of our spectra, photon- 
noise ( Bouchy et al.| 2001 ) contributes only marginally to the uncertainty of our Hermes RVs. 
Instead, the precision of our RVs is dominated by the intra-night stability of the wavelength 


^http://cds.u-strasbg.fr/ 

^see http;//obswww.unige.ch/$\sim$udry/std/std.html 
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BJD - 2400 000 

Phase 

Vr 

[ kms“^] 

a{vr) 

[ kms“^] 

55816.469627 

0.82074 

0.802 

0.018 

55816.470920 

0.82098 

0.813 

0.018 

55816.472221 

0.82122 

0.824 

0.018 

55817.510667 

0.01474 

-21.125 

0.018 

55817.511961 

0.01498 

-21.202 

0.018 

55818.494094 

0.19800 

-31.354 

0.018 

55818.495384 

0.19824 

-31.362 

0.018 

55818.496675 

0.19848 

-31.342 

0.018 

55819.497269 

0.38494 

-21.989 

0.018 

55819.498564 

0.38518 

-21.984 

0.018 

full table available at CDS 



Table 1: Sample of the new Hermes RV data that are published in their entirety at the CDS, 
and shown here for guidance regarding their form and content. Radial velocities have been 
shifted to the CORAVEL-ELODIE zero-point. Pulsation phase is dehned here as 0 when 
Vr = at the steep part of the RV curve. Uncertainties are hxed at 18ms“^ to account for 
uncertainty in the wavelength calibration and zero-point calibration, see text. 


HD 

Sp.Type 

(Vr) 

[kms“^] 

'^r,ref 

[kms“^] 

[kms“^] 

RMSne 

[ms“^] 

RMScorr 

[ms“^] 

10780 

KOV 

2.771 

2.70 

0.071 

38.7 

17.1 

32923 

G4V 

20.627 

20.50 

0.127 

70.1 

14.8 

42807 

G2V 

6.068 

6.00 

0.068 

115.2 

29.7 

82106 

K3V 

29.784 

29.75 

0.034 

45.8 

20.8 

144579 

G8V 

-59.452 

-59.45 

-0.002 

42.7 

20.0 

168009 

GIV 

-64.581 

-64.65 

0.069 

36.4 

20.0 

197076 

G5V 

-35.413 

-35.40 

-0.013 

78.3 

10.9 

221354 

KOV 

-25.111 

-25.20 

0.089 

53.8 

15.5 


Table 2: List of RV standard stars with spectral types from SIMBAD, Hermes RVs of stan¬ 
dard stars from new observations, Udry et al. (1999b), offsets between new and reference 
RVs, and RMS of new observations without (subscript nc) and with pressure corrections 
applied (subscript corr) following Anderson (2013, Sec. 2.1.5). 
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calibration and the long-term stability of the instrnment. 


To improve RV precision and stability, we employ a method developed in |Andersoii 


(2013, Sec. 2.1.5) to correct for the RV drift indnced by atmospheric pressnre variations that 
occur during the night. Pressure variations are the leading cause for drifts in Hermes RVs, 
since the temperature of the instrument is stabilized to within 0.01° K (Raskin et ah 2011). 
The RV drift correction due to pressure variations is based on changes in the refractive in¬ 
dex of air and has been shown to be precise to approximately 10ms“^ in the case of the 
Coralie spectrograph]^ for which RV drift corrections are measured using interlaced simul¬ 
taneous ThAr exposures. For Hermes, our method improves RV precision by a factor of 
approximately 2.5, as measured by the decrease in RMS for all standard stars. 


Figure]^ shows both the uncorrected (red dashed lines) and corrected RVs (blue solid 
line) for the standards stars HD 144579 and HD 168009 as a function of time for the different 
observing runs during which they were observed. The highest precision of 9m s“^ is achieved 
for HD 168009 during a ten-night observing run in July 2013. We find no evidence for long¬ 
term variations over the 3 year duration of the observations. Based on all standard star 
measurements, we estimate the long-term precision of pressure-corrected Hermes RVs to be 
approximately 15ms“^. Accounting for the additional uncertainty due to zero-point offsets 
(see above), we adopt 18ms“^ as our error budget for the investigations based purely on 
Hermes RVfl Note that this adopted uncertainty of 18 ms ^ is more than a factor 2000 
smaller than the pulsation-induced peak-to-peak RV amplitude of 38.6kms“^. 


3. Results from Spectroscopy 


In this section we describe our analysis of 6 Cephei’s RV curve, starting with the initial 


discovery of J Cephei’s nature as a spectroscopic binary (Sect. 3.1). After subtracting the 


RV drift due to orbital motion observed in Hermes data, we establish a high-precision ref¬ 
erence model for the pulsations alone, which we subtract from a combined data set of RVs 


from Hermes and several literature sources (cf. Sect. 3.2). Using this combined dataset, we 


determine the orbit for the binary system in Sect. 3.3 


We also searched the individual high-quality Hermes spectra as well as the cross¬ 
correlation prohles for a signature identifying the companion, hnding none. We note that 
previous analyses of lUE (Evans et ah 1993) and HST/COS ( [Engle et ah 2014) spectra did 


^Coralie is mounted to the Swiss 1.2m Euler telescope located at La Silla Observatory, Chile 
^However, the uncertainty adopted to determine JCephei’s orbit is 47ms“^, cf. Sect.and Fig.j^ 
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Fig. 1.— Hermes RV standard stars HD 144579 and HD 168009. The measurements are 
separated into individual observing runs and plotted as a function of time centered around 
the median observation date of a given star in a given observing run. We plot the mea¬ 
surements as residuals around the median of all measurements (i.e., including data from all 
observing runs) as listed in Tab.|^ We offset measurements from different observing runs 
and different stars (black dotted line) from one another for clarity. Pressure-corrected (see 
text) measurements are shown as solid blue connected dots, the uncorrected measurements 
are shown as a red dashed line. The labels on the right indicate observing run and RMS of 
the pressure-corrected RVs. 

































not report evidence of a companion’s signatnre, effectively rnling ont early-type compan¬ 
ions. Given the high signal-to-noise ratio (S/N) of onr spectra and the even better S/N 
of the cross-correlation prohles, we estimate that the companion must be at least a factor 
100 fainter than 6 Cephei, at least in optical bandpasses. Based on Geneva stellar evolution 
models (Georgy et al. 2013) and assuming an approximately 5Mq Gepheid, this implies a 
companion mass below approximately 1.75 Mq. This is consistent with the upper limit on 


companion spectral type (A3, i.e., M 2 < 2Mq) set by non-detection in lUE spectra (Evans 


1992). 


3.1. Discovery of <5 Cephei B using HERMES RVs 


Figurej^ shows the Hermes RVs obtained for 5 Gephei, phase-folded with the best-htting 
pulsation period of Ppuis = 5.366274 d. The ten-fold RV uncertainty is shown in the upper 
right corner. The observation date is traced by symbol color going from red (hrst observations 
in September 2011) to yellow (September 2014). The color coding reveals the average radial 
velocity to be time-dependent, thus revealing orbital motion caused by a hidden companion. 


We model the RV variations due to pulsation using the technique described in [Anderson 


et al. (2013). In a nutshell, we £t a Fourier series to model the phase-folded RV curve and 


increase the number of harmonics until an F-test indicates spurious £t improvement. This 
approach yields £t residuals that clearly exhibit a strong temporal correlation with an RMS 
of 337ms“^, see Fig.[^ 


By running the htting algorithm while assuming a model composed of the sum of a 
Fourier series and linear, quadratic, and cubic trends, we hnd that the orbital drift seen in 
Hermes RVs is best described by a cubic polynomial and that the best-£t pulsation model 
has 14 harmonics. We retain this model as our pulsation reference model for the following 
steps. 


Figure|^shows the residuals from Hermes RVs after subtracting the model that accounts 
for pulsations as well as the cubic drift due to orbital motion. The residuals are flat over 
the observational baseline. However, the phase-folded residuals do exhibit some structure, 
which is exposed by applying a color scale to trace the observation date. The RMS of our 
Hermes residuals is 47ms“^, i.e., higher than the 18 ms“^ estimated from RV standard stars 
in Sec. 12. II 


This higher-than-expected residual may be explained by small stochastic variations in 
the pulsation period seen in photometry of other Gepheids obtained with the Kepler and 
MOST satellites (so-called period-jitter, see Derekas et al. 2012 Evans et al. 2015), which 
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Fig. 2.— Phase-folded HERMES radial velocity curve with Ppuis = 5.366274 d. Observation 
date is traced by a color scale and increases from red to yellow. The ten-fold mean uncertainty 
of the measurements is shown in the top right corner. The RV offset at constant phase is 

due to the spectroscopic binary nature of S Cephei. 

1.0 

^ 0.5 

E 

^ 0.0 
u 

6 -0.5 
- 1.0 

55800 56000 56200 56400 56600 56800 57000 

BJD - 2400000 

Fig. 3.— Residuals for Hermes RVs minus the pulsation model as a function of observation 
date. The trend (dashed line) reveals the presence of 6 Cephei’s spectroscopic companion. 
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has been explained as being dne to snrface convection and grannlation (Neilson & Ignace 


2014). We note also that Anderson (2014) recently discovered RV cnrve modulation in 


Cepheids, albeit at much lower amplitudes. Since these effects (period-jitter or modulation) 
limit our ability to precisely reproduce the pulsation curve, we adopt the RMS value of 
47ms“^ as our Hermes RV uncertainty for 5 Cephei for the remainder of our analysis. 


3.2. Combination with literature data 


The cubic drift seen in the residuals in Fig.[^ indicates orbital motion at a timescale 
longer than the observational baseline achieved by our Hermes observations. We therefore 
searched the literature for data suitable for determining the orbit of 6 Cephei. 


6 Cephei is one of the most-studied variable stars and several authors have previ- 


ously published RV data for it, including 

Shane (1958), Barnes et al. 1 

1987), 

Wilson et al. 

(1989), 

Butler (1993) 

Bersier et ah (1994) 

, Gorynya et al. (|1996), Kiss 

(1998 

, Storm et al. 

(2004), 

Barnes et ah 

(2005). Historically, 5Cephei’s RV curve has been thought to be well- 


understood, which may in part explain why its spectroscopic binary nature has gone unno¬ 
ticed for so long. The main reason, however, is that 6 Cephei’s orbital signature has a small 
RV amplitude, high eccentricity, and long orbital period, and thus requires high-precision 
velocimetry over an observational baseline spanning at least two years. For comparison, the 
various RV datasets available in the literature have typical observational baselines on the 
order of one year and do not have sufficient precision to detect binarity during this timeframe. 


Zero-point offsets must be corrected for when combining RV data from different instru¬ 
ments and the literature. 6 Cephei represents a particularly difficult case, because zero-point 
differences can be on the same order of magnitude as the difference in due to orbital 
motion between datasets. Therefore, a well-determined common zero-point is the key to de¬ 
termining the orbit of 6 Cephei accurately. To this end, we adopted the CORAVEL-ELODIE 
RV zero-point (cf. Udry et ^|1999a|[b and Tab.|^. 

Of the available literature sources, we excluded the following from our analysis due to 
insufficient precision (cr^^ > 1kms“^): Shane (1958), Barnes et ah (1987), Wilson et ah 


(1989). We further discarded the measurements by Butler (1993) and Storm et ah (2004), 


since no information is available to determine the zero-point differences with CORAVEL. 


Finally, we excluded the single measurement published by Gorynya et ah (1996) 


Conversely, the measurements published by Bersier et ah (1994) and Barnes et ah (2005), 
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Fig. 4.— Hermes RV residuals (RMS= 47ms“^) after accounting for both pulsation and 
the cubic drift (see text) due to binarity. The left panel shows residuals as a function of the 
observation date, the right panel as a function of pulsation phase, with Ppuis = 5.366274 d. 
We trace observation date by color scaling the data points (red is oldest, yellow is newest) 
to expose additional signal that is not correctly modeled by pulsation and binarity and is 
likely related to random fluctuations in pulsation period (period-jitter). 


0-C [km/s] 
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as well as publicly available RV data from the ELODIE archiv^ ( [Moultaka et al. 2004) are 
already on the common CORAVEL-ELODIE RV zero-point. We found that an offset of 
—0.35kms“^ to the measurements by Kiss (1998) is appropriate to improve agreement with 


the (contemporaneous) data from the ELODIE archive and Barnes et ah (2005). This offset 
is similar to the precision stated by Kiss (1998, ~ 0.3kms“^). The determination of the 


zero-point offset between Hermes and CORAVEL-ELODIE is discussed in Sect. |3.1| above. 
Having thus calibrated the zero-point differences based on measurements of standard stars, 
we can use the combined data set to determine the orbit with conhdence. 

Classical Cepheids are known to exhibit changing periods due to their secular evolution. 
However, Cepheids can also exhibit erratic changes of unknown origin in their pulsation 


periods (see e.g. Berdnikov et al. 2000). When combining RV data from the literature, we 


noticed that variable periods have to be accounted for. We attempted to use the ephemerides 


and rate of (pulsation) period change by Berdnikov & Ignatova (2000) to obtain accurately 
phase-folded RV curves, but we were unable to obtain a satisfactory result. 

Since obtaining good phase-folding is required in order to correctly subtract the pul¬ 
sation reference model, we phase-folded the combined data set in the following way. First, 
we separated the dataset into three parts with different pulsation periods. The motivation 
for separating the evolution of the pulsation period in this way is the observation that the 


period of 5 Cephei changes very slowly (Eddington 1919, Berdnikov & Ignatova 2000). We 


then determined the best-£t pulsation periods for each of these three epochs by minimiz¬ 
ing the scatter in the residuals after fitting for the pulsation alone. We thus adopt the 


following pulsation periods: 5.3657 ± 0.0013 d for data by 

Bersier et al. 

(1994 

mean epoch 

JD 2444467.12); 5.36615 ± 0.0005 d for data from ELODIE, 

Kiss 

(1998 

, and 

Barnes et al. 


(2005, mean epoch JD 2450 398.61); 5.366274 ± 0.00006 d for Hermes RVs (mean epoch JD 
2456 430.38). While the center values of this sequence would imply a slowly increasing pul¬ 
sation period {dP/dt ~ 0.5 — 1.1 x 10“®syr“^), these adopted periods agree to within their 
uncertainties. After applying our updated pulsation periods, we shifted all three epochs for 
0 = 0 to occur at minimum radius, i.e., when the velocity is equal to on the steep part of 
the RV curve. This provides us with an accurately phase-folded pulsation curve from which 
we subtract the pulsation reference model to reveal the orbital motion of 6 Cephei. 

Figurej^ shows the RV curve based on the combined dataset. The residuals clearly 
demonstrate the presence of orbital motion. 


'http://atlas.obs-hp.fr/elodie/ 
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BJD - 2,400,000 

qP C? C? qP c? cP 

CN^ 

^ ^ ^ ^ <r^ 4^ 


Pulsation phase 
0.0 0.2 0.4 0.6 0.8 


1.0 



BJD - 2,400,000 


Fig. 5.— The combined dataset [Hermes + literature RVs) modeled only for pulsation- 
induced variability. The left panels show the measurements (top) and residuals (bottom) 
against observation date, the right hand panels show the same measurements against pulsa¬ 
tion phase. Observation date is traced by the applied color scale, with the newest measure¬ 
ments drawn white. 
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3.3. Orbit Determination using HERMES and Literature RVs 


Figurej^ shows the orbital motion of 6 Cephei exposed by subtracting our pulsation 
reference model from the phase-folded combined dataset. As can be seen in the top right 
panel, the orbit is well-sampled except for the ascending part of the orbital RV curve. Our 
observations did not sample this part of the RV curve, since it occurred shortly (a few 
months) before the start of our observations. 

Table[^ provides the orbital solution for 6 Cephei determined from the orbit-only com¬ 
bined RV curve using a standard Keplerian model. We used the tool Yorbit (Segransan, 
et al., in prep.) to hrst determine an initial orbital estimate via a genetic algorithm and 
then characterize the parameter uncertainties using the marginal distributions of Markov 
Chain Monte Carlo simulations with 500 000 iterations, cf. Fig.[^ From the orbital solution 
and assuming Mscep ~ 5.0 — 5.25 Mq (cf. Sec.[^ we determine the minimum mass of the 
companion to be 0.2 ± 0.02 Mq. While a^ez and Oi sin i are listed assuming Mscep = 5.25 Mq, 
the stated values remain within the stated uncertainties if Mscep = 5.0 Mq is adopted. 


Comparing our results to other known Cepheid orbits listed in the Cepheid binary 
database by Szabados (2003), we hnd that the semi-amplitude K of this orbit is the second 
smallest among all known binary Cepheid orbits, albeit with much larger eccentricity and 
somewhat longer period than W Sgr’s orbit (Groenewegen 2008). This also helps to explain 
why 6 Cephei has not previously been identihed as a spectroscopic binary. 


4. Searching for an Astrometric Orbital Signature 


Projecting the orbit derived in Sect. 3.3 to the distance determined by Benedict et al. 


(2007, w = 3.66 ±0.15mas) yields an orbital relative semimajor axis of 0^1 = 21.2mas, with 


Farm 

Unit 

Vj 

[kms“^] 

To 

[d] 

Povh 

[d] 

K 

[kms“^] 

e 

OJ 

[deg] 

®rel 

[au] 

oi sin i 
[10“^ au] 

0 

1 

CO 

±_, ^ 

Value 

-16.787 

55649.68 

2201.87 

1.509 

0.674 

246.77 

5.82 

226.3 

0.784 


0.026 

24.681 

5.73 

0.239 

0.038 

2.37 

0.18 

2 I. 9 I 

O. 249 I 

CT_ 

0.049 

19.861 

6.31 

0.080 

0.021 

4.90 

0.19 

7 . 9 I 

O.O 83 I 


Table 3: Orbital solution for 6 Cephei based on the combined Hermes and literature radial 
velocities. (T+ and (T_ denote the upper and lower standard errors derived from marginal 
distributions. Quantities with superscript dagger have been computed using Gaussian error 
propagation. Other uncertainties were estimated by the MGMG analysis. 
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Fig. 6.— Orbital solution for 5 Ceptiei based on the combined dataset from which we sub¬ 
tracted our pulsation reference model. We trace observation date by color scaling data points 
from red for the oldest to yellow for the newest measurements. Left panels show measure¬ 
ments (top) and residuals (bottom) against observation date, right panels show the same 
data against orbital phase, assuming Porb = 2201.87d, see Tab.[^ and 0 = 0 at pericenter 
passage. We indicate the range of orbital phases at which Hipparcos, HST, and Gaia have 
observed or will observe 6 Cephei. 
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an orbital barycentric semimajor axis of the Cepheid of Oi > 0.84mas, a value more than 


hve times the uncertainty estimated by 

Benedict et ah 

(2(1(17 

), and even seven times the 

uncertainty stated by 

van Leeuwen et ah 

(2007 

, Hipparcos, 3.71 ± 0.12mas). This prompts 


the question whether the observations taken by Hipparcos or HST are sensitive to the orbital 
motion, cf. Fig.j^for information as to which ranges of orbital phase were observed by these 
missions. We therefore explore the sensitivity of Hipparcos and Gaia astrometry to S Cephei’s 
hidden companion in this section. 


4.1. Hipparcos Intermediate Astrometric Data 


To test whether Hipparcos (Perryman & ESA 1997) was sensitive to the orbital motion 
of S Cephei, we analyze the Hipparcos intermediate astrometric data (lAD) published on the 
DVD attached to the new reduction by van Leeuwen (2007[). 


4-1.1. Parallax from Intermediate Astrometry Data 

A 5-parameter £t to the 95 measurements given on the DVD yields the residuals shown in 
Fig.[^ left panel. There are several outliers and £t quality is poor (x^d = 2.01, RMS(0—C) = 
1.43mas). We obtain a best £t parallax of = 4.37 ± 0.27mas. 

We therefore discarded a total of 6 data points on the basis of their excess residual from 
the following satellite orbits: 180, 252, 396, 759, 1786, and 2126. We then repeated the 5- 
parameter £t to 89 measurements and obtained a better £t (Fig.[^ right panel; = 0.65, 
RMS(0—C) = 0.81 mas), as well as smaller parallax: rtJsparm = 4.09 ± 0.16 mas. Note that 



Fig. 7.— Marginal probability density functions for orbital period, eccentricity, and semi- 
amplitude as determined by Yorbit’s MCMC algorithm. The dashed line indicates the 
median value adopted as numerical reference. 
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there is a discrepancy of AciJ = 0.38 mas (0.21 mag) between the parallax valne published in 


the van Leeuwen (2007, w = 3.77 ± 0.16 mas) reduction and our result based on the lAD. 


There may be several reasons for such a discrepancy, including: 


1. Our htting routine. We excluded this possibility by processing other stars, in particular 
5 Cephei’s visual companion HD 213307 for which we obtained a parallax of 3.69 ± 0.46 
mas in exact numerical agreement with the published value (3.69 ± 0.46 mas. 


van 


Leeuwen 2007). As an additional cross-check, we applied the same methodology to 


more than 20 other Cepheids using the same routine and obtained results that agree 


well with the published values of van Leeuwen (2007). 


2. our selection of lAD. As noted in the lAD header, 3% of the lAD were discarded to 


obtain the solution of van Leeuwen (2007), although it is not specihed which ones. We 


attempted to reproduce their result by discarding only 3 measurements but did not 
succeed. 


3. the cluster solution. It appears that van Leeuwen (2007) used a special procedure 


for cluster stars to better remove outlier measurements for individual stars, under the 
assumption that all cluster stars are at indistinguishable distance. While our parallax 
result is larger than the previously published value, it is well within the wide range of 
parallaxes (2.57 — 5.28 mas) reported for presumed members of CepOB6. 

4. other unknown procedures or modihcations that may have been applied specihcally to 


van Leeuwen 

(2007 

) and 

van 


Leeuwen et ah (2007). 


Our parallax estimate {w = 4.09 ± 0.16 mas) is also considerably {2a) larger than the 
HS'T-based result by Benedict et ah (2002, w = 3.66 ± 0.15mas). Such an increase in 


parallax would increase 5 Cephei’s absolute magnitude by 0.24 mag, making it intrinsically 
fainter than previously thought {My = —3.23 mag using the absolute magnitude published 
by Benedict et al.|2002). We note that the iLS'T-based period-luminosity relations presented 


by Benedict et ah (2007) seem to agree better with such an increase in absolute magnitude 


for 5 Cephei. 

The HST observations were taken at an orbital phase at which a signihcant parallax bias 
due to orbital motion is not very likely, cf. Fig.|^ However, the assumption of a physical 
association between HD 213307 as well as 6 Cephei in the loose association Cep OB6 was 


required to “reduce [the HST] astrometric residuals to near-typical levels” (Benedict et ah 


2002 , Sec. 5.5). 
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Concerning 6 Cephei’s membership in Cep OB6, we note that the association’s discovery 
did not take into account radial velocity data, since only little such information was avail¬ 
able at the time (de Zeeuw et ah 1999). Inspection of the SIMBAD database yields radial 
velocity information for 10 of the 19 presumed member stars (in addition to 6 Cephei), with 
values ranging from —7kms“^ to —38kms“^ and median uncertainty of 1.2kms“^ based 


Fehrenbach et al.f ( 

1996) 

Grenier et al. 

(1999) 

Famaey et al. 

(2(1(15 


Gontcharov (2006); Kharchenko et ah (2007). The RVs are distributed as follows: HIP110807 
and HIP112998 have ~ —7±2kms“^, HIP 110497 has = —13.3±0.5kms“\ HIP109492 
and HIP110988 (= HD213307) are close to 6 Cephei’s = —16.8kms“^(to within lkms“^), 
HIP113993 has Vr = —20.9 ± l.lkrns"^, HIP113316 has Vr = —25.6 ± 8kms“^, and three 
stars (HIP110266, HIP110275, HIP110356) have Vr < —30kms“^ with reported uncertain¬ 
ties between 3 and 15kms“^. This relatively wide range of RVs suggests that CepOB6 is 
not gravitionally bouncQ i.e. that the assumption of a common distance for all presumed 
member stars is not valid. This interpretation is corroborated by the wide range of parallax 
values of the presumed member stars (see item 3 above). However, more homogeneous and 
high-quality radial velocity measurements of the presumed members of Cep OB6 are required 
to further illuminate this issue. 


Although the above suggests that not all of Cep OB6’s presumed members have indis¬ 
tinguishable distance, we note that the observational evidence does support the physical 
association between HD 213307 and 6 Cephei, since their parallaxes and radial velocities 
agree to within the uncertainties. A re-assessment of the HST astrometric data without 
the assumption of cluster membership (this affects e.g. the spectrophotometric parallax for 
HD 213307, which is part of the HST reference frame), and taking into account 5 Cephei’s 
orbital motion would be useful for testing whether the difference between the published HST 
parallax and our result can be reconciled. 


4.1.2. Orbit Analysis 


We search for an orbital signature in Hipparcos lAD using the methodology described in 


Sahlmann et al. 

(2011b 

), which has been shown to reliably detect such orbital signatures in 

Hipparcos lAD ( 

Sahlmann et al.||2011b|a 

Sahlmann & Fekel|2013 

). We use the spectroscopic 


orbital parameters given in Tab.to £t the lAD with a seven-parameter model, which has 
the free parameters inclination i, longitude of the ascending node D, parallax vu, and offsets 


®For this to be the case, the RV dispersion of member stars should not exceed a few kms ^ (e.g. 


Mathieu 


19861. 
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to the coordinates (Aa, A5) as well as offsets to the proper motions (A/ia, A/i5). We 
then search a two-dimensional grid in i and hi for its global x^-minimnm with a nonlinear 
minimization procedure. We determine the statistical significance of the derived astrometric 
orbit via a permutation test ( Zucker fc Mazeli||2001 ) for which we employ 1000 pseudo-orbits 
and derive parameter uncertainties using Monte Carlo simulations that include propagation 
of RV parameter uncertainties. 


Only 41% of the orbit are probed by Hipparcos measurements, whose average measure¬ 
ment uncertainty of 1.07 mas is furthermore greater than the barycentric minimum semima¬ 
jor axis of the orbit (0.84mas). It may therefore not surprise that we determine virtually 
identical results for 5-parameter (single star) and 7-parameter (binary) models, with flat 
residuals even for the 5-parameter model (Fig.j^. From an F-test, we obtain a probability 
of 4.4% that the single-star model is true, and the permutation test yields an orbit detection 
signihcant at the l.Scr level (93.7%). While the evidence for the orbital signature is only 
marginally significant, it has been determined with two independent methods that are in 
agreement. It thus appears that there is some temporal coherence present in the Hipparcos 
data, although we cannot claim orbit detection using Hipparcos astrometry. We can, however, 
use our orbital analysis to set (very) loose constraints on the inclination (10° ^ i 170°) 
of the orbit, see Fig.|^ allowing us to place an upper limit on the companion mass with 
0.2 Mq < M 2 < 1.2 Mq, which is fully consistent with the lack of spectral features due to 
the companion in the Hermes data set, cf. Sect.|^ 

While we are unable to claim detection of the orbit from Hipparcos astrometric data, 
we caution that the previous estimate of 6 Cephei’s proper motion may have been affected 
by the companion. 


4.2. Gaia 


The ESA space mission is currently conducting an unprecedented census of our 

Galaxy, measuring position, proper motion, and parallax for more than a billion stars during 
a nominal mission duration of five years. 6 Cephei will be among these billion objects thanks 
to GaiaA ability to observe very bright starG*^ (Martm-Fleitas et ah 2014). For an assumed 


single-measurement precision of croaia ~ 100 /ras, the minimum barycentric semimajor axis 
(oisinz 840/ras (Sect. 3.3)) is roughly 8 times larger than the measurement uncertainty. 


"http://www.cosmos.esa.int/web/gaia 

Stars brighter than G-band magnitude 5.7 will be heavily saturated and need special treatment to 
determine their centroids. We conservatively assume a single observation per Gaia field-of-view transit. 
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Figure]^ shows the range of orbital phases covered during the nominal mission duration plus 
a possible one year extension. It shows that Gaia measurements will cover the majority 
of the orbit and, crucially, will be measuring astrometry during periastron passage. We 
investigate the detectability of 6 Cephei’s orbit from Gaia astrometry in two ways. 


First, the detectability of astrometric orbits with Gaia can be estimated from a consid¬ 
eration of the astrometric signal-to-noise S/N = ai y/N /acaia (Sahlmann et al.[2011a, 2015). 


Assuming 86 Gaia observations of 5 Cephei during the 5-year mission (see below) and ac¬ 
counting for 10% dead time with an assumed 100 p.as accuracy for individual measurements 
and the minimum amplitude of ~840 /ras, we obtain S/N ^ 73, which is much higher than 


the detection threshold of 20 described in Sahlmann et al. (2015). This indicates that Gaia 


will detect the astrometric orbit of 6 Cephei, despite not measuring at all orbital phases. 


Second, we use the Gaia astrometric simulation software AGISlab (Holl et al. 2012) 


to investigate the detectability of <5 Cephei’s orbit derived in Sect. 3.3 We adopt the Ax 


metric from Perryman et al. (2014) to evaluate the detectability of the astrometric binary 


signal. This metric measures the reduction in minimum x^ when going from a (single-star) 
5-parameter solution to a (binary) 12-parameter Keplerian solution. Even for our adopted 


worst-case scenaric^^, we hnd an improvement of AxminSo- = 139 when accounting for the 
orbital motion. This is well above the threshold for precise (parameters determined to 
better than 10%) orbit characterization (Perryman et al. 2014, Ay^ > 100). We therefore 


predict that Gaia will clearly detect and characterize the astrometric counterpart to the 
spectroscopic orbit discovered here. 


5. Discussion: Piecing Together the Puzzle 


Now that 6 Cephei’s nature as a spectroscopic binary is revealed, it is worth revisiting 
other observed features of the prototype of classical Cepheids in this new light. 


Engle et al. (2014) recently provided evidence that d Cephei is a soft X-ray source with a 
luminosity of Lx(0.3 —2keV) ~ 4.5 — 13 x 10^® ergs“^ and peak flux at kT = 0.6 —0.9 keV. As 
these authors discuss, young (~ 120 Myr for a 5 Mq Cepheid, according to Geneva evolution 
models by Ekstrom et al.||2012; Georgy et al.||2013) low-mass main-sequence companions can 


provide coronal X-ray emission. In Sects. pA] and |4T721 we constrained the mass range for 
the unseen companion to be 0.2 < M 2 < 1.2 Mq based on radial velocities and Hipparcos 


^^This assumes the worst possible configuration of sky-alignment and observation noise, as well as 20% 
dead time 
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astrometric measurements. As Cepheids are rarely detected in X-rays, it appears likely that 
6 Cephei’s young main sequence companion is responsible for the detected variable X-ray 
activity. 

Due to the high eccentricity of the orbit, 6 Cephei and its companion have recurrent close 
encounters, with a pericenter distance of rper = (1—e)arei = 1.89 au = 409 Rq = 9.5 i?*, where 
i?* = 43.3 dTurn erl 


1988 Merand et ah 2005; Natale et ah 2008). It is also important 


to bear in mind that 5 Cephei is currently in the core He burning phase, most likely on 
the second crossing of the instability strip as shown by the decrease in pulsation period 


(measured e.g. by E. Hertzsprung and reported by Eddington| 1919, Berdnikov & Ignatova 


2000, and Engle et ah 2014). It has therefore previously occupied the red giant branch, 
where its radius was even larger, approximately R-^^rg ~ 80 i?©. Assuming an unchanged 
orbit, pericenter passage would have brought the two stars to within 5 x 

Taking these considerations further, we consider the possibility of previous interactions 
between 6 Cephei and its companion. Since 6 Cephei is on a highly eccentric orbit and the 
Roche-lobe formalism is valid only for circular orbits, we adopt the Roche-lobe formalism in 


the quasi-static approximation as presented by Sepinsky et ah (2007a). Assuming Mscep = 
5.2 Mq, M 2 = 0.7Mq, and equatorial velocity Veq = 10kms“^, we obtain the volume- 


2.8 R 


*,RG- 


equivalent Roche lobe radius at pericenter of Rnoche.corr = 1.05 au ~ 5.2 R* 

Hence, at pericenter passage S Cephei hlls 19% of its (quasi-static) Roche lobe, and the 
Roche lobe radius at pericenter is approximately 55% of the distance at pericenter passage 
(^"per = 1.89au). At apocenter, the distance between the two stars is a factor (1 -|- e)/(l — 
e) = 5.1 larger, and 5 Cephei hlls only about 3.7% of its Roche lobe. 5 Cephei may thus 
become noticeably deformed due to tidal interactions close to pericenter passage, while being 
spherical near apocenter passage. During the red giant phase, this situation would have 
been even more extreme, if the past orbit was similar to the present-day orbit. The above 


considerations suggest that the observed circumstellar material (Marengo et ah 2010 Merand 
et al.||2006 ) and bow-shock (Matthews et ah 2012) may originate from previous and ongoing 
binary interactions. A more detailed investigation of such interactions is required to discuss 


this scenario in terms of 5 Cephei’s mass loss history (see the discussion in Matthews et ah 


2012 ), but is considered out of scope for this paper. 


If 6 Cephei and its companion have a history of episodal interactions at pericenter pas 
sage, then the high eccentricity of the orbit may appear surprising. However, Sepinsky 


et ah (2007b) showed that rapid circularization is not expected for all close-in binaries with 
this type of interactions. Applying our results for 6 Cephei to their formalism shows that 
6 Cephei’s orbit is not expected to have been rapidly circularized. Furthermore, 6 Cephei is 
a visual binary whose outer companion is understood to be physically associated and on a 
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very long-period orbit (Benedict et al. 2002). The high eccentricity of the inner binary (the 


one shown in the present work) conld thus have been driven up by the outer companion by 


a Kozai-Lidov mechanism (Lidov 1962; Kozai 1962) and may have varied significantly over 
its evolutionary history. 

If 6 Cephei and its companion have undergone significant interactions, one might expect 
the evolutionary status of 6 Cephei to vary significantly from that of a Cepheid with a single 
star progenitor. We therefore examine 6 Cephei’s current evolutionary status in Figure 10 


to search for signs of non-standard evolution. We compare observed absolute V magnitude, 
{B — V)o color (both from Benedict et al.|2002), and rate of (pulsation) period change Ppuis = 


of Solar metallicity that include rotation (lEkstrom et al. 

2012 

Georgy et al.||2013D 

and have 

been studied specifically in the context of classical Cepheids by 

Anderson et al. 

(2014 

). Strong 


disagreement between the observations and these predictions could be considered evidence for 
binary interactions, since no binary interactions are accounted for in these model predictions. 
However, we find that 5 Cephei’s location in both diagrams is consistent with a 5.25 Mq 
C epheid whose progenitor had a slightly faster-than-average initial rotation (fl/flcrit ~ 0.7). 
Assuming these parameters for the evolutionary models yields an age of 112 Myr. Adopting 
our new parallax from Sec. 4.1.1| yields a smaller initial mass of M ^ 5.0 Mq and thus an older 
age of 127 Myr, while leaving the implications regarding the rotational history unchanged. 


At the present level of accuracy, we thus do not hnd any irreconcilable discrepancies 
between the predicted and observed evolutionary states of 6 Cephei. This shows that bi¬ 
nary interactions, if present, have either had a negligible effect on the evolutionary path 
of 6 Cephei, or that the interactions are weak and slow enough for 6 Cephei to reach an 
equilibrium state similar to a non-interacting star. 


In summary, the discovery of 6 Cephei’s nature as a spectroscopic binary helps to com¬ 
plete the puzzle created by a plethora of observations. While there is no clear evidence 
for a non-standard evolutionary path, 6 Cephei is a particularly interesting example of the 
limitations that this evolutionary phase is subject to for binary stars (cf. 
and deserves detailed observational follow up and dynamical modeling to further investigate 
its intriguing past that may have been marked by tidal interactions due to both the inner 
(discovered here) and outer (HD 213307) companions. 


Neilson et al.||2014) 


^^Models accessible at http;//obswww.unige.ch/Recherche/evoldb/index/Interpolation/ 
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6. Conclusions 


230 years after the discovery of its variability (Goodricke 1786), we discover the spec¬ 


troscopic binary nature of 5 Cephei, archetype of classical Cepheid variable stars and one of 
the most-studied variable stars. 

Our discovery is demonstrated using new high-precision radial velocities measured from 
high-quality optical spectra obtained with the high-resolution spectrograph Hermes. Com¬ 
bining these new high-precision data with lower-precision RVs from the literature, we deter¬ 
mine the orbital solution for the spectroscopic binary, which is an inner binary to the outer 


visual binary system discussed by Benedict et ah (2002). 5 Cephei thus appears to be a pair 
of binary stars. 


We re-analyze Hipparcos intermediate astrometric data (cf. Sects. 4.1.1) and obtain the 


parallax tu = 4.09 ±0.16 mas {d = 244 ±10 pc) using a 5-parameter (single star) model. This 


result is larger than the estimates reported by van Leeuwen (2007, w = 3.77±0.16 mas). 


van 


Leeuwen et al. 

(2007 

w = 3.71±0.12 mas), and 

Benedict et al.|( 

2002 


these previously published results based on Hipparcos and HST astrometry agree to within 
their stated uncertainties, they all shared a common assumption of S Cephei’s membership 
in the loose association Cep OB6, which our analysis does not and which we argue should be 
revisited using radial velocities. Relaxing this assumption and accounting for orbital motion 
in a re-analysis of the HST astrometric data would be useful to test whether the existing 
HST astrometry can be reconciled with our Hipparcos-heised result. 


We perform an orbital analysis of the Hipparcos lAD in Sect. 4.1.2 and hud tentative 


evidence for an orbital signature, although no detection can be claimed. Based on detailed 
simulations, we show that Gaia is highly sensitive to the astrometric orbit of 5 Cephei and 
will likely model the full set of Keplerian parameters with better than 10% accuracy. The 
orbit will have to be accounted for when determining proper motions from Hipparcos and 
Gaia data. 

Using the constraints provided by the optical spectra, the orbit measured from RVs, the 
astrometric orbital analysis, and assuming a mass of 5.0 — 5.25 Mq, we constrain the mass 
range of the companion to be 0.2 < M 2 < 1.2 Mq. Adopting the lower mass for 5 Cephei 
mainly affects the upper mass limit, which would become 1.1 Mq in this case. Given that 
the spectroscopic companion is expected to be the same age as S Cephei, i.e., approximately 
100—130 Myr (depending on mass and ZAMS rotation rate of the progenitor), the reported X- 
ray emission detected using XMM-Newton (Engle et al.|2014) could be explained by magneto- 


rotational activity of a young main-sequence star. 


The close periastron approach of the two stars has potentially far-reaching consequences 
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for the explanation of the observed circnmstellar environment of 6 Cephei. Detailed modeling 
of the orbital and stellar evolntion of this complex system is desirable to farther improve onr 
nnderstanding of the archetype of classical Cepheids and its intrigning past. 
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Fig. 8.— Hipparcos fit residuals for all 95 data points (left panel) and after removing 6 
outliers (right panel). Time is relative to epoch J1 991.25 as in the original Hipparcos catalog 
(Perryman & ESA 1997| . 
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Fig. 9.— Joint confidence contours on the i —fJ-grid. Contour lines correspond to confidence 
levels at la (solid), 2a (dashed), 3a (dotted), and 4a (dash-dotted). Crosses indicate the 
position of the best non-linear adjustment solution for each of the 100 Monte Carlo samples 
of spectroscopic parameters; the star corresponds to the adopted parameters. 
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Fig. 10.— 6 Cephei’s color-magnitude and color-rate of period change diagrams that compare 
the measured quantities with predictions from Solar metallicity Geneva stellar evolution 
models ( [Ekstrom et al. 2012 Georgy et ah 2013; Anderson et ah 2014). Note that the 
predicted rates of period change are based on the evolution of the average density, i.e., they 
have not been determined from a pulsation code. We use absolute magnitude and dereddened 
color by Benedict et al.| (2002) and the rate of period change measured by Engle et al. (2014) 
as the observed values. The left panel is particularly well-suited for determining progenitor 
mass, whereas the rate of period change is very sensitive to the main sequence rotation of 
6 Gephei’s progenitor. The legend in the right panel indicates ZAMS mass in solar units 
followed by the model’s initial rotation rate, where VI corresponds to a; = G/Gcrit = 0.1, 
V7 to o; = 0.7, and V9 to a; = 0.9. V5 corresponds to w = 0.568, i.e., v/vcnt = 0.4, see 


Ekstrom et al. (2012). Gomparing the observed values to model predictions (assume single 


star evolution), 5 Gephei is consistent with a progenitor mass of approximately 5.25 Mq and 
slightly faster-than average surface rotation [u ~ 0.7). Adopting our new parallax estimate 
(tu = 4.09 ± 0.16) would change My to —3.23 mag, resulting in a lower inferred mass of 
5.0 Mq while yielding the same result in terms of rotation 

































